IS A CURVED FLIGHT PATH IN SAR BETTER THAN A STRAIGHT ONE? 
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Abstract. In the plane, we study the transform -R 7 / of integrating a unknown function / over 
circles centered at a given curve 7. This is a simplified model of SAR, when the radar is not directed 
but has other applications, like thermoacoustic tomography, for example. We study the problem 
of recovering the wave front set WF(/). If the visible singularities of / hit 7 once, we show that 
the "artifacts" cannot be resolved. If 7 is a closed curve, we show that this is still true. On the 
other hand, if / is known a priori to have singularities in a compact set, then we show that one can 
recover WF(/), and moreover, this can be done in a simple explicit way, using backpropagation for 
the wave equation. 



1. Introduction 

In Synthetic Aperture Radar (SAR) imaging a plane flies along a curve in R 3 and collects data 
from the surface, that we consider flat in this paper. A simplified model of this is to project the 
curve on the plane, call it 7; then the data are integrals of a unknown density function on the 
surface over circles with various radii centered at the curve. Then the model is the inversion of the 
circular transform 



(1) R 1 f( r ,p)= f(x)d£(x), p€ 7 , r>0, 

J \x— p\=r 

where d£(x) is the Euclidean arc-length measure, and the center p is restricted to a given curve 7(i). 
This transform has been studied extensively; injectivity sets for i? 7 on have been described in 
full [3], see also [7]. In particular, each non-flat curve, does not matter how small, is enough for 
uniqueness. In view of the direct relation to the wave equation, this transform, and its 3 dimensional 
analog, see section |4j have been studies extensively as well and in particular in thermoacoustic 
tomography with constant acoustic speed, see, e.g., [H [21 [TOj, [ITJ [T2l [21 [151 US]- A related 
transform is studied in [H [6] . 

The problem we study is the following: what part of the wave front set WF(/) can we recover? 
Clearly, we can only hope to recover the visible singularities: those conormal to the circles involved 



in the transform, see also section 3.1 



If 7 is a straight line, there is obvious non-uniqueness due to symmetry. Moreover, we can 
have cancellation of singularities symmetric about that line. More precisely, we can recover the 
singularities of the even part of / and cannot recover those of the odd part. 

Based on this example, it has been suggested that a curved trajectory 7 might be a batter flight 
path. This question has been studied in [17J, and some numerical examples have been presented 
suggesting that when the curvature of 7 is non-zero, the artifacts are "weaker", and with increase 
of the curvature, they become even weaker. By artifacts, they mean singularities in the wave front 
set of R*R~ff that are not in WF(/) located at mirror points, see Figure [2j The same problem but 
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formulated in terms of the wave equation model problem has been studied from a point of view 
of FIOs in [18] , see also [8] , where the artifacts have been explained in terms of the Lagrangian of 
R~. They found that the artifacts are of the same strength, as an order of the corresponding FIO. 
More precisely, this is true at least away from the set of measure zero consisting of the points whose 
projections to the base falls on 7 (points right below the plane's path, i.e., r = 0), and for (x, £) such 
that the line trough it is tangent to 7 at some point. The latter set is responsible for existence of a 
submanifold of the Lagrangian near which the left and right projections are not diffeomorphisms. 
What part of the singularities of / can be recovered however has not been studied, except for the 
cases when there is an amplitude which vanishes at the mirror points; then the artifacts can be 
ruled out by a prior knowledge. 

The main purpose of this paper is two fold. First, we study the local problem — what can be said 
about WF(/) knowing WF(i? 7 ) near some point, which localizes possible singularities of / near 
two mirror points. More generally, we can assume that each line through WF(/) crosses 7 once, 
transversely. Then we show in Theorem 2.1 that curved trajectories 7 are no better than straight 
lines — singularities can still cancel; moreover, the artifacts are unitary images of the original. We 
also describe microlocally the kernel of i? 7 modulo C°°. For simplicity, we stay away from the 
measure zero set mentioned above. While this can be generalized globally for arbitrary curves, 
without the single intersection condition, we do not do this but study a closed curve encompassing 
a strictly convex domain. Then we show again that recovery of singularities is not possible. In this 
sense, a curved or even a closed path is no better than a straight one. 

On the other hand, when 7 is closed and strictly convex, if we know a priori that WF(/) lies 
over a compact set (i.e., the projection of WF(/) onto the x-space is in a fixed compact set), then 
we show in Theorem |2.3| that one can recover the visible singularities. We even present a simple 
way to do that in the interior of 7, by backprojecting boundary data for the wave equation, see 



Proposition 3.3 In this sense, a curved trajectory is better. The effect which makes it possible 
is based on the fact that any singularity inside should be canceled by two outside if we see no 
singularities on the boundary; but the latter should be canceled by other singularities farther away, 
etc. At some point, this sequence would leave the compact set over which WF(/) lies a priori, thus 
contradicting the assumption on /. 

This transform belongs to the class of the X ray transforms with conjugate point studied by the 
authors in [21]. The circle centered at 7 and passing through x in the direction has a conjugate 
point at the mirror image of (x, #~ L ). The approach which we follow here is different however. 
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2. Main Results 

Fix a smooth non self-intersecting curve (si, S2) B s *— > j(s). For convenience, assume that s is an 
arc-length parameter. We parameterize FL,f then by s and the radius r > 0, so we write i? 7 /(r, s) 
instead of i? 7 /(r, 7(5)), compare with 0. The possible obstruction to recovery of singularities is 
well understood. Fix an orientation along 7 by choosing the normal field 7 := (—72,71)- This 
defines a "Left" and a "Right" side of 7 near 7. Let (xl,£l) £ T*R 2 \ 0, and xl 7- Assume 
that the line through (xl,£l) intersects 7 form the left, at some point po = 7(so)> and that this 
intersection is transversal. If it is tangent, then the Lagrangian of i? 7 is not of a graph type, see 
e.g. [18]. We call such a singularity visible from 7. We want to emphasize now that visible does 
not necessarily mean recoverable from i? 7 /, which is the whole point of this paper. Let xr be 
the point symmetric to xl about the line tangent to 7 at p (a "mirror" point w.r.t. p), and let 
£r be the symmetric image of £7,, see Figure [2] Note that £l, ^r may both point towards 7, or 
both point away from it. Then (xl,£,l) and (xr,£,r) are symmetric images to each other w.r.t. the 
symmetry about that tangent line lifted to the cotangent bundle. Denote this symmetry map by 
C, i.e., C(x l ,±Cl) = (xr,±£r). 

Set to = \xl~ 7(so)| = \xr—j(so)\. The circular transform i? 7 /(r, s), for (t, s) close to (to, sq) and 
acting on a function / supported in a small neighborhood of xl and m can only detect singularities 
close to ±£l and ±£r respectively, see section |3.1[ but it is not clear if it can distinguish between 
them. We can expect that a singularity at (xl, might be cancelled by a singularity at (xr, ±£r) 
and we might not be able to resolve the visible singularities. 




Figure 2. Mirror points: (xr,±£r) are mirror points to (:e£,±£l), and vice versa 



Any open conic set in T*R 2 \ satisfying the assumptions so far (also implied by the assumption 
below), can be written naturally as the union S^U Tir of two sets satisfying 

For any E Sl (or S_r), the line through (x,£) hits 7 

transversely from the left (right), at exactly one point different from x. 

Condition Q implies that Ex, Hr are unions of disjoint open sets: Ex = E^ U E^, E^ = E^ U E^ 
where the positive and the negative signs indicate that x+t£ hits 7 for t > and t < 0, respectively. 
Then C : E^ — > E^. Let / be a compactly supported distribution with WF(/) C E^ U E#. The 
question we study is: what can we say about WF(/), knowing WF(i2 7 /)? Since i? 7 is linear, it is 
enough to answer the following question: let R~f G C°°(7 x R+) (or let be smooth microlocally 
only, in a certain conic set). What can we say about WF(/)? 
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Without loss of generality, we can assume that C(El) = 'Er. In section 3.4 below, we show that 
R-y, restricted to distributions with wave front sets in Hl,r is an FIO associated with a canonical 
graph denoted by Cl,r, respectively. In particular, the projection it(Cl(xl, on the base is 

(to, so), i.e., to is the time it takes to get to 7 with unit speed, and sq corresponds to the point po 
where that line hits 7. Then we set 

(3) S 7 :=C L (S L ) =C R (Z R ). 

The possible singularities of Rjf with / as above can only be in E~. 

Theorem 2.1. Assume Q). Let Yjl, Hr, £ 7 , C be as above, and let fi, Jr be compactly supported 
distributions with WF(/i) C Ex, WF(/r) C T,r. Then there exists a unitary Fourier Integral 
Operator U with canonical relation C so that 

(4) R^fL + fn) GC°°(S 7 ) <=> fR-Uf L €C°°p R ). 
Moreover, U = —AJ^Al, with Al.r described in section 3.4 and Proposition 



3.2 



The unitarity of U above is considered in microlocal sense: U*U — Id an UU* — Id are smoothing 
in Si and T>r, respectively, where the adjoint is taken in L 2 sense. 

The practical implications of Theorem 2.1 is that under assumption ([2]), only the singularities 



of f R — Ufi (or, equivalently, U*Jr — fi) can be recovered. We can think of it as the "even part" 
of / in this case. In particular, for any fi £ V' with WF(/i) C E/, there exists Jr E V with 
WF(/h) C Y,r so that -R 7 (/l + Jr) G C°°(E 7 ). An explicit radial example illustrating this is 
presented in Example [l| Thus the artifacts when using R*Ryf to recover WF(/) are not just a 
problem with that particular method; they are unavoidable, and they are a unitary image of the 
original, i.e., "equal" in strength. From that point of view, a curved path is no better than a 
straight one. 

We study next the case where the path 7 and / are such that there are singularities (x,£) of / 
for which the line through them hits 7 more than once. Of course, this can happen for a curved 
path only. Consider the examples in Figure [3| where each of the dashed lines intersects 7 at most 
twice. We assume that there are no more intersection points than shown. On the left, the trace 
that (x,£) leave on 7 at p\ can be canceled by its mirror image (37, £1) about (the tangent at) p\. 
Equivalently, (27, £1) can create an artifact at (27, £1), and vise-versa; related by a unitary map. 
Similarly, the singularity on 7 caused by (x, £) at p_i can be canceled by its mirror image (x-\, £-1) 
about We assume there that the lines through (27, £1) and do not intersect 7 again. 

If we know that one of the three singularities cannot exists, then none does. In particular, we can 
recover (x,£) if we know a priori that either (27, £1), or (x_i,^_i) cannot be in WF(/). Without 
any prior knowledge, we cannot. On the right, all those five singularities can cancel if they are 
related by suitable unitary operators. If we know that one of them cannot be in WF(/), then none 
can. 

Notice that P-1PP1P2P3 is a geometric optics ray reflected by 7. To obtain (2^, £2), for example, 
we start from going along the broken path, and at any point between P2 and P3, we go back 

from the same distance but along a straight line. If we go along the broken ray past p^ (not shown 
on the picture), and come back along a line the same distance, we end up at X3. The point x~\ 
can be obtained similarly, going in direction opposite to £. 

So far we assumed that each line appearing in the construction intersects 7 at most twice. If 
this is not true, the mirror points to (x, £) form a directed graph. We will not study this case. 

Assume now that 7 is a closed curve and it encompasses a strictly convex domain fi. The 
discussion above suggests the following. For any (x,£) G T*R 2 \ 0, let <&*(x,£) = (x + i£/|£l>£)- 
Let be defined on T*fi \ in the same way for small \t\, then extended by reflection, etc. At the 
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values of t corresponding to reflections, we take the limit from the left. We call this path, extended 
for all positive and negative t, a broken line through (x,£). Then all mirror points of (x,^), where 
possible artifacts might lie, are given by 

(5) M(x,Z) = {$- t o$fa,£); teK}, xefl. 

This is a discrete set under our assumption, see, e.g., |16j . In the examples in Figure [3j this set is 
finite in each case, consisting of (x,£), (x±i,£±i), etc. Since 7 is a closed curve now, in our case it 
is infinite. The next theorem says that if we have a priori knowledge that would allow us to rule 
out at least one of those artifacts, then we can recover a singularity at Otherwise — we 

cannot. 



A 



A 






Figure 3. Singularities that cannot be resolved. Left: has mirror images 

and Singularities at any two of those three points are related by 

unitary maps. Right: an example with more than three points. 

Those arguments lead to the following "propagation of singularities theorem" . 

Theorem 2.2. Let 7 = dQ, where C R 2 is a strictly convex domain. Let f G T>'(R?) with 
7Hsupp/ = 0, and assume that Ryf G C°°. Then for any (x , f) G T*tl\0, either M(x,£) C WF(/) 
or M{x, n WF(/) = 0. 

As in the example above, if we know a priori that one of those points cannot be in WF(/), then 
none is, and in particular, / is smooth at One such case is when WF(/) a priori lies over a 

fixed compact set. 

Theorem 2.3. Let 7 be as in Theorem^^ Let f G £'(R 2 ). If R y f G C°° , then WF(/) contains 
no singularities visible from 7, and in particular, /|q G C°°. Moreover, f\^ can be obtained from 
Rryf modulo C°° by the back-projection operator described in Proposition 3.3, 



If we do not have a priori information about /, then WF(/) cannot be reconstructed. 

Then there is f G £>'(R 2 \ 7) \ C° 



Theorem 2.4. Let 7 be as in Theorem 



2.2 



so that Ryf G 

C 00 (R+ x 7). Moreover, for any f with singsupp/ C 0,, there is g with singsuppg C R 2 \Q so 
thatR^(f-g)£C co (R + x-/)/ 

The second statement of the theorem says that we can take any / singular in Q, and extend 
it outside so that its circular transform will be smooth on 7. Therefore, not only singularities 
cannot be detected but any chosen in advance / singular in £1 can be neutralized by choosing 
suitable extension singular outside O. We refer also to Example [T] and the remark at the end of it 
for a radial example. 

Those problems and the methods are related to the thermoacoustic problem with sources inside 



and outside O, see Remark 3.1 . 
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3. Proofs 

3.1. The wave front set of the kernel of R~. The Schwartz kernel of i? 7 is given by 

K(r,s,x) = -S s i 

r \ r 

The factor 1/r is not singular for r > 0, where we work. By the calculus of wave front sets, 

WF(^) = |(r jS ,x,dF* Sja: r ? ); x- 7 ( S )erS\ 7? = ^"^ , k ^ oj , 
where F = (x — 7(s))/r. Set uj = F\ s i to write this as 

WF(ft) = l(r,s,j(s)+ru,-^,-^u-j(s),*Lj) ; w G S l , r > 0, A; ^ o| . 
Set £ = —kco/r; then u; = — = e l£l> where e = ±1 is the sign of k, to get 

(6) WF(ft) = | (r, s, x, -e|e|, e ■ 7(«), ~0 ! x + ^ y| = 7{s),r > 0, e = ill . 

By the calculus of wave front sets, if we invert the sign of the sixth component there, £, and consider 
WF(7£) as a relation, this tells us where WF(/) is mapped under the action of Rj. Comparing this 
with the definition (18), (19) of Cl, and similarly for Cr below, we get 

WF( J R 7 /) c C L (WF(/)) uC fl (WF(/)) 

for / such that for any (x,£) G WF(/), the line x + s£ through (x, £) meets 7 exactly once, for 
s / 0. 

Let (r, a) be the dual variables to (r, s). The reason we use r instead of the more intuitive choice 
p for a dual variable to r is that by applying the 'I'DO A below, we will transform r into a variable 
denoted by t. Then 

(7) WF(iV)c{M<|T|}. 
Moreover, by ([2]), for any such /, we have |<r| < S\t\, S < 1. 

3.2. Reduction to a problem for the wave equation. Let u solve the problem 

{dj - A)u = inRjX R|, 

(8) < u|t=o = 0, 

d t u\ t =o = f, 

and set Af = m|r +X 7, i.e., 

(9) Af = Si ^P/ • 

\L>\ R-+X7 

The well known solution formula then implies 



rHyf(r, s) 



(10) A/(t,a) = J J^LJ dr, p G 7. 

Our assumptions imply that Rjf(r,p) = for < r < 1. The integral above then it has a kernel 
singular at the diagonal t = r only. It belongs to the class of Abel operators 

(11) Ah(t)= f4=dr, t>0. 
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Then 

(12) A = (A®Id)R r , 

in other words, A is just ARj but A acts in the first variable. The explicit left inverse of A is 

(13) h(r) = Bh(r) :=—■$- f . h ® dt, r > 0. 

Proposition 3.1. The operator A restricted to £' (R+) is an elliptic ^DO of order —1/2 with 
principal symbol 

a p (B)(r,r) = y/^e'^^ (r+ 1/2 + irl 1/2 ) . 

77ie operator B on £'(R+) is an elliptic ^fDO of order 1/2 principal symbol given by the 
inverse of that of A. 

Proof. The Schwartz kernel of A can be written as 

A(t,r) = A$(t,r,t - r), where A*{t,r,w) 



r 



y/t + r 

and 104. = max(ro, 0). The Fourier transform of w^ 1 ^ 2 is equal to 

/- -in/ 4 ( -I/ 2 , ■ ~V 2 
/7re ' [ Ti + ir_ 



-1/2 



Then A is a formal ^DO with an amplitude given by the partial Fourier transform of A^ w.r.t. w, 
i.e., 

W4 r / -1/2 . -1/2 



V* + 



Since t and r are strictly positive, there is no singularity in 1/y/t + r. The singularity at ^ = can 
be cut off at the expense of a smoothing term. Set t = r to get the principal symbol of A. Since B 
is a parametrix of A, the second assertion follows directly. □ 

Note that the full symbol of A can be computed from the asymptotic expansion of the Bessel 
function Jo since A is the composition of the Fourier Sine transform and the zeroth order Hankel 
transform, see 



3.3. Working with the Darboux equation. The unrestricted spherical means Gf(t,x) := 
(27rt) _1 Rf solve the Darboux equation 

(df + ^d t -A\ Gf(t,x) = 

with boundary conditions Gf(0,x) = f(x), dtGf(0,x) = 0, see e.g., [2] and the references there. 
The Darboux equation has the same principal symbol as the wave equation and therefore the same 
propagation of singularities for t ^ 0. Replacing the wave equation with the Darboux one seems as 
a natural thing to do — this would have eliminated the need for the operators A and B. On the 
other hand, t = is a singular point which poses technical problems with the backprojection, and 
for this reason we prefer to work with the wave equation. 
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3.4. Geometric Optics. The solution of ([8]) is given by 

sin(t|£>|), e-' ll \ D \ e«l D l 

(14) u= — , , f = : — r /H : — -f = U + + U-, 

K ' \D\ 2i\D\ 2i\D\ J + ' 

where 

_ 1 f j((x-v)-e-t\i\) 



(15) 

1 ' ,i((*-»)-€+*lil)_L 



■/(y)dyd^. 



- (2vr)2y 2i|e|' 

The first term u_ is in the kernel of dt + and if we consider t as a parameter, it is an FIO 
associated with the canonical relation (x, £) H >• (x + i£/|£|,£)- The second term u- is in the kernel 
of dt — i\D\ associated with (x, £) h-» — £)■ 

We assume now that WF(/) C S L , see Then we set A L f = Af with / as above. We define 
Ar in a similar way. 

Restrict Q to R x 7, see @, to get 

(16) A L = A+ + A L , 

where A±f are the restrictions of the two terms above to R x 7. For the first term, we set x = j(s) 
to get 



(17) A+f := - e ——f „ = (2vr)- 2 / eKWHJ^tDi^/^djide 



-it|B| 



2i|D 



-2 / P i((7(*)-s/)-c-*iei) 



(-1) 



R+x 7 v ' y 2iiei 



Since we made an assumption guaranteeing that 7(5) • £ 7^ 0, this is an elliptic FIO with a non- 
degenerate phase function, see e.g, [22], Ch.VI.4 and Ch.VIII.6, of order —1 associated with the 
canonical relation 

ct ■ (7(«) - te/ieu) — ► (*, -lei, m • o. 

well defined on £j. Another way to write this is the following. Let t(x,£) > 0, s(x,£) be such that 
a? + t(x, = 700>£))- Then 

(18) C+ : (x, — >• (t(x, 0,s(x, 0, -1^, 7(s) " 0- 

Similarly, the second term in ( |14[ ) defines 

1<|D| (2vr)- 2 |e i ((^)^)^l)^L/( y )d y de 



R+x 7 

This is an FIO associated with the canonical relation 

(19) C L : (x, — >• (t(x, -0, s(x, -0, |e|,7(s) ■ £), 

since x — t(x, — £)£/|£| = j(s(x, — £)) for € ££. We now define as on Sj, and on 

S L . Similarly, A L / is defined as A+/ when WF(/) G £+ Also, set £ = C(£ L ) C T*(R+ x 7), see 
also @. 

We define C^, Cr, A^/, A/j/ in the same way. In fact, they are the same maps as the "L" ones 
but restricted to £^, £#, and / with wave front sets there, respectively. Clearly, the map C defined 
in the Introduction satisfies 

C = C R Cl ■ Si — ► £r, 

and ^ holds. 

Relations (18), (19) imply also the following, compare with 0, 

(20) WF(A± „/) C {(t, s, r, a); \a\ < T Sr} , 
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where < 5 < 1 is the cosine of the smallest angle at which a line through (x, £) G WF(/) can hit 



7, see (18), (19) 



Since and Ar are elliptic FIOs (associated with canonical graphs), they have left and a right 
parametrices A^ 1 and A^ 1 , of order 1 associated with C^ 1 and C^ 1 , respectively. We have the 
following more conventional representation of those inverses. 

We recall the definition of incoming and outgoing solutions in a domain f2. Let u(t, x) solve 
the wave equation in [0, T] x £2 up to smooth error, i.e., {df — A)u G C°°, where H C R 



2 



IS 



a fixed domain, and T > 0. We call u outgoing if u(0, ■) = 1^(0, •) = in f2; and we call u 
incoming if u(T, •) = ut(T, •) = in f2. We micro- localize those definitions as follows. A solution 
of the wave equation modulo smooth functions near R x 7, on the left (or right) of 7 is called 
outgoing/incoming, if all singularities starting from points on R x 7 propagate to the future only 
(t > 0), and respectively to the past (t < 0). 

Proposition 3.2. Let ul be the incoming solution of the wave equation with Dirichlet data h on 
R_i_ x 7, where WF(h) C E; and assume Q). Then 

(21) A L 1 h = 2d t u L \t= . 

Proof. Call the operator on the r.h.s. of (21) M for a moment. To compute MAif, recall 



Assume first that WF(/) C Sj. Then Aif = Aj/, see (17), i.e., Aif is the trace on the boundary 
of -u+ defined in ( |15[ ). Now, to obtain MA^f, we have to find the incoming solution of the wave 
equation with boundary data A^f. That solution would be u+ modulo C°°, i.e., u+ = ul in this 
case. Then MA^f = 2dtu+\t=o, by the definition of M. The latter equals /, by the definition of 
u + . If WF(/) C S^, then Aif = A^f, and MAif = 2dtU-\t=o = f ■ In the general case, / is a 
sum of two terms with wave front sets in and E7, respectively. 

To see that M is a right inverse as well (which in principle follows from the characterization of 
Ax as an elliptic FIO of graph type), let ul be as in the proposition. Then Mh = 2dtUi\t=Q. To 
compute A^Mh, we need to find first the outgoing solution of the wave equation with Cauchy data 
(0, Mh) at t = 0. This solution must be ul- Indeed, call that solution v for the moment and write 



v = v+ + V- as in ( 14). Assume first that WF(/i) is included in r < 0, where r us the dual variable 
to t, see(|18[). Then the singularities of v+ hit R x 7 but those of V- do not, by ([2]). The solution 
v + has Cauchy data at t = given by 

(22) {~\D\~ x Mh^Mh\ 



see (|14j). Now, ul has the same Cauchy data, which proves that u+ = ul- Then A^Mh is the 
trace of u+ on the boundary, which is h. The case r > 0, and the general one, can be handled in a 
similar way. □ 

3.5. Proof of Theorem [2TT} Set f = f L + fn- Assume now that Ryf G C°°(£). Apply A ® Id 
to that, where A acts w.r.t. to r and Id is w.r.t. s, to get Af G C°°(S). Since ^4 has a left inverse 
on R + , this is actually equivalent to R^f G C°°(S), i.e., 

(23) ^/eCT) ^ Af L + Af R GC°°^). 



Indeed, recall that (r, a) is the dual variable to (t, r); then ^4 £B) Id is elliptic on {r / 0}. By (20), 
A Id is elliptic in a conic neighborhood of WF(A/), which proves our claim. The restrictions of 
the wave front sets of fi, and f R imply that we can replace A above by its microlocalized versions 
Al )jR : 

(24) Ver(S) Ai/i + Afi/fiGrfE). 
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Now, apply the parametrix A^ 1 to (24) to get 
(25) /L+A^erp). 



Of course, starting from (25) we can always go back to (|24j). Therefore, (13) and (24) are equivalent, 
and they are both equivalent to A^AlJl + Jr £ C°°(S). 

To show that U is unitary, we will compute H^/lII first. Let fi be as above. Denote by ul the 
solution with Cauchy data (0, fi) at t = 0. To obtain A^ 1 , we need to solve backwards (to find the 
incoming solution) of the wave equation on the right of 7 with boundary data A+/ = U£|r, + X7 . Let 
us call that solution ur. On the other hand, ul restricted to the right of 7 is an outgoing solution 
with the same trace on the boundary. Then v := ur — ul solves the wave equation the right of 
7 , and for t = we have v = ur; while for t = T > 1, we have v = —ul- Moreover, v has zero 
Dirichlet data on the boundary. Therefore, up to a smoothing operator applied to /l, the energy 
of ul at t = T coincides with that of ur at t = 0. The former one is equal to the energy of the 
Cauchy data (0,/l), up to smoothing operator, and therefore, E(ur(0)) = + K /l||i,2, where 
K is smoothing. If WF(_/x) C SJ, then ul solves {dt + \\D\)ul £ C°°, and then so does ur. Then 
E(ur(0)) = HA^j 1 ^!!^, see (22), where h = Al/l- Therefore we showed that 



\\Q.d + K)f L \\v = ||A^A L / L || i2 . 

This proves that U*U = Id modulo an operator that is smoothing on Sj. In the same way we 
show that this holds on which is disconnected from Since U is microlocally invertible on 
Si, we get that U is unitary up to a smoothing operator on E^, as claimed. 
This completes the proof of the theorem. 

Example 1. We give an example of cancellation of singularities. Let 7 = S 1 be the unit circle 
parameterized by its polar angle s. Then I7I is not unit but constant, which is enough. Let / 
be the characteristic function of the circle |x| = 1/2, i.e., f(x) = H(l/A — \x\ 2 ), where H is the 
Heaviside function. Then, clearly, Ryf(r,9) is singular at r = 1/2 (not only), with a singularity of 
the type y/r — 1/2, see also Figure We will construct a radial function g supported outside the 
unit disc so that Rj(f — g)(r, 0) is smooth in a neighborhood of r = 1/2. 

We will work with radial functions only, i.e., functions of the form Fd^l 2 ). We will identify the 
latter with F, somewhat incorrectly. Then R^F is independent of the angle s, and it is enough to 
fix s = corresponding to x = (1,0). Then we have 



l R~ f F(r) = I F ((1 + r cos 6) 2 + (r sin 6) 2 ) d9. 

T J —TT 



The factor 1/r can be explained by the requirement that the measure along each circle is Euclidean. 
Since 1/r is a smooth factor near r = 1, we will drop it. We also use the fact that the integrand is 
an even function of 9, so we denote (l/(2r))i? 7 by R: 



rn 

Rf(r)= f (l + r 2 + 2rcos9) d9. 
Jo 



Set r = 1/2 + h. We are interested in the singularities near h = and in what follows, \h\ 1. 
After replacing 9 by ir — 9, we get 



p7T 

Rf(l/2 + h)= H(2rcos9 - r 2 -3/4) d# = H(h) arccos 

70 r=l/2+h 



(1/2 + hf + 3/4 
l + 2/i 
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The following calculations were performed with Maple. The series expansion of the expression 
above is 

(26) Rf (1/2 + h) = H(h) (V2h^ - ^V 2 + + 0(h^ . 
We are looking for a radial g of the type 

(27) g(\x\ 2 ) = H(t)(a + a x t + a 2 t 2 + ... )| Hx |a_ 9/4 . 
Then 

rn 

Rg{\/2 + h)= H(h)(a + a x t + a 2 t 2 + . . ■)\ t={l+2h)cose+{l/2+h) 2^ / ^0 

J 

= a A (h) + axAxQi) + a 2 A 2 (h) H . 

For Aq we easily get 



By (26), to cancel the /i 1 / 2 term in R(f — g), we need to chose 

a = V3/3. 



Then 

~6 

To improve the smoothness near h = 0, we compute 



" 5o) = -^r^l 2 + 0(h 7 / 2 ), g := a H(t)\ t=lxl 2_ 9/4 . 



Al(h) = ^3/2 _ 37^2^5/2 + 0{h 7/2y 

3 15 



Then, as before, we find that we need to choose 

a 2 = -5/16 

to kill the 0(h 3 ^ 2 ) term, and then 

83^ 
T20 



R(f~ 9i) = -^7T^ 5/2 + 0(h 7 / 2 ), g x := H(t)(a + ai t)\ t=lx]2 _ 9/4 . 



Note that this was possible to do because the leading coefficient (the one in front of /1 3 / 2 ) in the 
expansion of A± is non-zero. The latter also follows from the ellipticity of A. Proceeding in the 
same way, we can get a full expansion of the conormal singularity of g at | x | = 3/2 that would 
make R(f — g) smooth at r = 1/2. 

The first three coefficients of g are shown below 

We could continue this process to kill all the singularities for all r, not just at r = 3/2 by con- 



structing a suitable jump of g at r = 5/2, then at r = 7/2, etc., which also illustrates Theorem 2.4 
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Figure 4. Left: density plot (white = 1, black = 0); right: a graph of / and g with 
R(f — g) smooth near r = 1/2. 
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Figure 5. The thick line: The graph of R(f—g) near r = 1/2 computed numerically 
with three terms in the expansion of g. The dotted line: the graph of Rf having a 
square root type of singularity. 



3.6. Proof of Theorem |2.2[ Let first 6 WF(/). Let x( x >D) be a pseudo-differential 

cutoff in some small conic neighborhood of (x, £). Then R^x( x iD)f has singularities in small 
neighborhoods of C^x, £) in T*(R + x 7), where we used the notation above. Take the plus sign 



first. Since R^f is smooth, there must be another singularity that cancels this one. By section 3.1 
it must be at the mirror point (27, £x) of (x,£) about the line tangent to 7 at pi, see Figure [2 
Clearly, (37, £1) belongs to the set M.(x, £); it corresponds actually to the first point t, when t, 
increases from t = 0, not equal to (x, £), see §5§. Moreover, by Theorem 2.1 we can construct f\ 
with a wave front set near that mirror point so that Rj(xf + /1) is smooth at C + (x,£). Since the 
line through (27, £1) crosses 7 transversely, it has to cross it again, also transversely. This creates 
another singularity on the boundary, represented in Figure [2] by P2 (of course, that singularity is an 
element of T*(R + x 7)). It needs to be canceled by another one, etc. We repeat the same argument 
for C~(x, £). Therefore, we showed that if (x,£) £ WF(/), then the whole set M.(x, £) is in WF(/). 

Now, assume that (x,£) WF(/). Since R-yf is smooth, a singularity at (37, £1) cannot exist, 
because it can only be canceled by one at (x,£). Then we show step by step that no point in 
M(x,£) can be singular for /. 



SAR 
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3.7. Proof of Theorem 2.3, Let (xcbCo) be visible from 7. This means that xo G" 7, and that the 
line through (xo,£o) intersects 7 transversely. Then (xq,£o) G M(x,£) for some (x,£) G T*f2 \ 0, 
even if xo is outside f2 (if xo G ^2, then (x,£) = (xoi £()))• Indeed, for this, we have to show that 
the equation o 3>*(x,£) = (xo,£o)> x G $7, is solvable for some i. This is true because we 
can set (x, £) = 3>^* o <3?*(xo, £o)> where t is any interior point of the non-empty interval (ii,t2) 
for which xo + i^o G SI. For any t, the projection 7r o $*(x, £) onto the base is in fi. Then 
1 7T o <1>~' o <3?^(x, £)| > \t\ — Ck, where Ck '■= max(|x|; x G Cl). Therefore, M(x,£) does not lie over 
any compact set. By the compactness assumption of the theorem, it has elements outside WF(/). 



Then by Theorem 2.2 (x, £) WF(/). 



3.8. Constructing a parametrix for R^f in Q,, when WF(/) lies over a fixed compact set. 

We will give another, constructive proof of Theorem |2.3| for the singularities of / inside T*£l. Let 
singsupp / C K , where K is a fixed compact set. Fix T so that 

(28) T > max(|x- y\; x G d£l,y G K) . 

Then all singularities of the solution u of ^ would leave Cl for t > T, and i? 7 / G C 00 for r > T. 
The latter is obvious even without the propagation of singularities theory. Let v be the incoming 
solution of the wave equation in fi with Dirichlet data A/ = (A Id)i? 7 / on [0, T] x (9f2, cut-off 
smoothly near t = T. More precisely, let x be a smooth function of t so that %(t) = for t > T, 
and = 1 for < t < To, where To < T is chosen so that Tq satisfies (28) as well. Let v solve 



(29) 



(dl - A)v 

u\t=T 

d t u\ t =T 
u \[o,T]xdn 



in [0, T] x O, 

in Q, 

in n, 
h. 



where h will be chosen in a moment to be x^-f ■ Set 



(30) 

Then G X ^f 
(31) 



/ in modulo C 

{ {d 2 t - A)w 

w\t=T 

d t w\t=T 
w \[o,T\xdQ 



Gh = d t v\ t =o. 
Indeed, consider w : 



u 



G 
G 





(1 - xK 



v. It solves 
in [0, T] x O, 



Then / — G/i = c?titf|t=o G C°°(Q), which proves our claim. 
To summarize this, we proved the following. 



Proposition 3.3. Let 7 be as in Theorem 2.2 and let f G V(R 2 ) be such that singsupp / C K \ r y, 
where K is a fixed compact set. Let T > 0, x be as in (2tfy and flsO. Then 



G X {A®Id)R 1 f = f\ Q modC°°(0). 

To complete the proof we only need to notice that by assumption, singsupp / is at positive 
distance to 7 = d£l, which guarantees that WF(/i), with h = x(A <X> Id)i? 7 /, is separated from 
t = 0, and the singularities of w are never tangent to dtt. This makes the operator G an FIO of 
order with a canonical relation a graph, like in the previous sections, and in particular G is well 
defined on such h. Therefore, Gx(A ® Id) R^f is well defined. 
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3.9. Proof of Theorem 2.4 We first present a proof along the lines of the proof of Theorem 2.2 
above. We prove a somewhat weaker version first: for any T > 0, we can complete / to a distribution 
in R 2 \ 7 so that R 7 f G C°°((0,T) x 7). Fix (s,£) G T*fi \ 0, and let / has a wave front set in 
some small neighborhood of that point. Let uq be the solution of the wave equation in the plane 



with Cauchy data (0,/). Then by section 3.1, R~f will only have singularities on T*7 near points 
defined by the line through (x, £) which lie over and (t\,pi), where t±i are the arrival 

times, see Figure |3j To cancel them, we chose g\ with singularities near (a;_i,£_i) and {x\,^\) 



see Figure [3] again, unitarily related to the singularity of / near see Theorem 2.1 Then 



+ u lj where u\ is the solution with Cauchy data (0, —gi), will have no singularities near the 
points mentioned above which project to (t-i,p-\) and (ti,p\). On the other hand, u\ will cause 
new singularities at points above (£±2>P±2)j see Figure [3] where only P2 is shown. We then construct 
g2 and a related 112 that would cancel them, etc. After a finite number of steps, the time component 
of the points {t±k,p±k) above which we have a singularity, will exceed T, and then we stop' and set 

g = gi + 52 + Then we use a microlocal partition of unity to construct g so that f — g would 

have the required properties without the assumption on WF(/). 

To prove the general case (i.e., to take T = 00 above), let (the subscript k now has a different 
meaning) be the distribution corresponding to T = n. Then g^ — g m = (f — g m ) — f — gp. has 
a circular transform smooth on (0, min(fc, m)) x 7, and gt — g m = in Q. The only possible 
singularities of that distribution could be those with the property that the line through each one of 
them intersects 7 transversely; then that singularity will leave a trace on 7. This implies that there 
are no singularities with travel time to 7 less than mm(k, m). Therefore, on some ball centered at 
the origin of radius min(/c, m) — C, the distribution g^ coincides with g m up to a smooth function. 
Then we can easily construct g as a "limit" of gk with a partition of unity, and this g would have 
the property i? 7 (/ - g) G C°° . 

Remark 3.1. The main results in this paper are also related to the thermoacoustic/photoacoustic 
model with sources inside and outside Q. The wave equation then is the underlying model and 



there is no need of the operator A. Theorem 2.1 and Theorem 2.4 then prove non-uniqueness of 



recovery of WF(/) as singularities of the data, with partial or full measurements. Theorem 2.3 



proves that this is actually possible if singsupp / is contained in a fixed compact set. The recovery 



is given by time reversal with T as in (28). The only formal difference is that in TAT, the wave 



equation is solved with Cauchy data (/, 0) at t = instead of (0, /); and the time reversal operator, 



see (21) and (30) does not contain dt- 



4. The 3D case: Recovery of the singularities from integrals over spheres 

CENTERED ON A SURFACE. 

Let r be a given smooth (relatively open) surface in R 3 . Let 
(32) Rrf(r,p) = f /(*) dS x , r > 0,p G T, 

J \x— p\=r 

where dS x is the Euclidean surface measure on the sphere \x — p\ = r. We show below that the 
results of the previous section generalize easily to this case as well. 

We assume again that / G £'(R 3 ) is supported away from T, and that for any (#,£) G WF(/), 
the line through (x,£) hits T once only, transversely. The main notions in section [3] are defined in 
the same way with a few minor and obvious modifications. In ^ and in the definitions (18), ( |19| ) 



of Cl,r we need to replace 7 • £; by the projection of £ onto the boundary, i.e., onto T*T, where 
p G r is the point where the line through (x, £) hits T. 



SAR 
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In this case, Rrf is more directly related to the solution of the wave equation; indeed 

u(t,x) = —R v f(t,x) 

is the solution of the wave equation in the whole space with Cauchy data (0, /) at t = restricted 
to R + x T. Then A = (4vrt)- 1 i? r , compare with §n§. Multiplication by (Ant)' 1 is, of course, an 
elliptic ^DO for t ^ (which is implied by our assumptions), and we get that Theorem 2.1 applies 
to this case, as well. In particular, we get that microlocally, we cannot distinguish between sources 
inside and outside the domain occupied by the 'patient's body" in thermoacoustic tomography. 
If the external sources have compactly supported perturbations, then we can, and time reversal 
reconstruct the singularities for a large enough time T such that each singularity coming from 
outside would exit before time t. This has been observed numerically in [15J. 

Finally, we remark that in applications to thermoacoustic tomography, the wave equation point 
of view is the natural one, actually. Then those results extend to variable speeds using the analysis 
in EDI. 
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